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Abstract 

During  the  sediment  acoustics  experiment  (SAX99),  a  number  of  independent  sound  speed 
and  attenuation  measurements  were  made  in  a  well-characterized  sandy  sediment.  These 
measurements  covered  a  broad  frequency  range  and  made  it  possible  to  test  both  Biot  theory 
and  Buckingham’s  more  recent  grain-to-grain  attenuation  model.  While  Biot  theory  was 
able  to  model  the  sound  speed  well,  it  was  unable  to  predict  the  attenuation  measured  above 
50kHz.  This  paper  presents  a  series  of  measurements  made  in  the  laboratory  on  a  simple 
glass  bead  sediment  to  test  the  hypothesis  that  this  deviation  was  due  to  scattering  within 
the  sediment.  The  sediments  used  were  saturated  with  two  fluids  with  different  viscosities 
in  order  to  shift  the  dispersion  into  the  frequency  range  of  the  measurement  system.  The 
measured  attenuation  in  the  glass  bead  sediments  exhibited  the  same  frequency  dependence 
as  observed  in  the  SAX99  data  indicating  that  scattering  is  not  the  relevant  attenuation 
mechanism.  A  sediment  model  is  proposed  which  is  a  hybrid  of  Biot  and  Buckingham’s 
theories.  This  hybrid  model  is  able  to  model  both  the  attenuation  and  sound  speed  in  the 
SAX99  sediment  and  in  the  laboratory  sediments.  The  implications  of  the  success  of  this 
hybrid  model  are  also  discussed. 

PACS  numbers:  43.30. Ma,  43.20. Jr 
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I.  INTRODUCTION 


In  the  study  of  sediment  acoustics,  there  is  an  ongoing  debate  about  the  nature  of  dis¬ 
persion  and  attenuation  in  unconsolidated,  sandy  sediments.  Many  measurements  have 
been  made  in  sandy  sediments,  both  in  situ  and  in  the  laboratory,  as  well  as  in  simpler 
unconsolidated  glass  bead  sediments.  The  results  of  these  experiments  have  led  to  two 
differing  conclusions.  Some  measurements  seem  to  indicate  that  attenuation  in  these  sedi¬ 
ments  varies  as  the  first  power  of  the  frequency,  /*,  with  little  or  no  significant  dispersion. 
This  view  has  found  it  s  most  vocal  support  in  Hamilton*  and  has  recently  been  revisited 
by  Buckingham  who  argues  that  this  conclusion  is  supported  by  current  research  into  the 
properties  of  dry  sand2  4.  Both  Hamilton  and  Buckingham  have  argued  that  the  sediment 
is  best  described  by  a  viscoelastic  model.  Both  of  these  approaches  determine  the  rele¬ 
vant  parameters  through  empirical  fits  to  the  data  limiting  their  predictive  value5.  The 
approaches  differ,  however,  in  that  Hamilton’s  approach  is  purely  empirical  while  Bucking¬ 
ham  believes  that  the  attenuation  is  due  to  friction  between  the  sand  grains4.  Buckingham 
feels  that  the  relevant  parameters  can  be  determined  from  the  fluid  and  sand  properties 
although  it  is  not  clear  at  this  point  in  the  theory’s  development  how  this  would  be  done. 

Due  to  the  porous  nature  of  the  sediment,  some  researchers  feel  that  the  relative  motion 
of  the  pore  fluid  and  the  sediment  frame  should  lead  to  viscous  damping  of  the  acoustic 
wave.  Biot  theory  provides  a  successful  description  of  sound  propagation  through  porous 
materials6’7  and  has  thus  been  applied  to  modeling  sand  sediments.  The  attenuation  pre¬ 
dicted  by  Biot  theory,  however,  varies  as  / 2  at  low  frequencies  and  /V2  at  high  frequencies 
which  contradict  some  of  the  measurements  in  sandy  sediments.  The  sound  speed  is  also 
predicted  to  exhibit  significant  dispersion  in  the  frequency  range  where  the  transition  be¬ 
tween  these  two  limiting  attenuations  takes  place.  Some  researchers,  most  notably  Stoll8, 
note  that  there  are  deviations  from  the  /*  frequency  dependence  in  many  measurements  as 
well  as  evidence  of  dispersion  indicating  that  viscous  losses  due  to  the  motion  of  the  pore 
fluid  may  be  important.  Indeed,  many  measurements  have  been  successfully  modeled  with 
Biot  theory  by  assuming  complex  frame  moduli  to  account  for  friction  and  local  viscous  flow 
at  the  grain  contacts.  This  modified  theory  has  come  to  be  known  as  Biot- Stoll  theory9. 
Biot  theory  has  an  advantage  over  the  viscoelastic  models  of  Hamilton  and  Buckingham 
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in  that  the  relevant  parameters  in  the  theory  can  in  principle  be  measured,  although  not 
without  some  effort.  Biot-Stoll  theory  remains  empirically  driven,  however,  in  that  the 
complex  frame  moduli  are  generally  determined  from  a  best  fit  to  the  data. 

With  the  number  of  measurements  available,  it  might  seem  strange  that  a  consensus 
has  not  been  reached  on  the  proper  description  of  sound  propagation  in  sediments.  The 
problem  lies  not  with  a  lack  of  experimental  effort  but  rather  in  the  difficulties  inherent 
in  performing  these  measurements.  Usually  it  not  possible  to  make  sound  speed  and  at¬ 
tenuation  measurements  over  a  frequency  range  broad  enough  to  observe  the  dispersion 
predicted  by  Biot-Stoll  theory  nor  is  it  usually  possible  to  fully  characterize  the  sediment. 
Realizing  the  importance  of  overcoming  these  difficulties,  a  number  of  researchers,  utilizing 
different  measurement  techniques,  came  together  to  determine  both  the  sound  speed  and 
attenuation  in  the  sediment  at  a  single  site  during  a  recent  high-frequency  sediment  acous¬ 
tics  experiment  (SAX99)10,11.  Extensive  efforts  were  made  to  determine  the  properties  of 
the  sediment  and  several  researchers  measured  the  sound  speed  in  the  sediment  over  a 
frequency  range  of  125Hz-400kHz  and  attenuation  over  a  frequency  range  of  2.6-400kHz. 
Details  of  these  measurements  can  be  found  in  Ref.  12. 

Despite  the  broad  range  of  frequencies  and  careful  measurement  of  the  sediment  proper¬ 
ties,  the  results  of  the  SAX99  experiments  have  not  been  successful  in  resolving  the  sediment 
acoustics  debate.  Both  Biot-Stoll  theory  and  Buckingham’s  model  have  been  applied  to  the 
data  and  both  have  their  respective  strengths  and  weaknesses  as  shown  in  Fig.  1.  Biot-Stoll 
theory  captures  the  dispersion  in  the  sound  speed  across  most  of  the  frequencies  with  the 
exception  of  the  two  lowest  frequency  measurements.  The  measured  attenuation  deviates, 
however,  from  that  predicted  by  Biot-Stoll  theory  at  the  highest  frequencies  and  appears 
to  follow  the  / 1  frequency  dependence  observed  previously.  Consequently,  Buckingham’s 
model  fits  the  attenuation  well  at  these  high  frequencies  with  only  a  slight  deviation  from 
the  lower  frequency  measurements.  His  theory  is  not  as  successful  in  modeling  the  sound 
speed  as  it  predicts  a  curvature  in  the  dispersion  that  is  opposite  that  of  the  measured 
values  and  only  fits  the  data  over  a  limited  frequency  range. 

As  Williams  et.  al.12  points  out,  these  results  are  one  of  four  separate  measurements 
that  indicate  that  the  sound  speed  in  sand  sediments  at  low  frequencies  is  lower  than  at 
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high  frequencies.  This  is  consistent  with  the  predictions  of  Biot  and  the  authors  speculate 
that  the  deviation  in  attenuation  at  high  frequencies  may  be  due  to  scattering  within  the 
sediment  from  either  volume  inhomogeneities  or  shell  and  shell  fragments.  In  order  to 
test  this  hypothesis,  a  series  of  laboratory  measurements  were  made  in  simple  glass  bead 
sediments.  The  measurements  were  performed  using  one  of  the  systems  which  was  used 
during  SAX99  to  measure  sound  speed  and  attenuation  over  the  100-260  kHz  frequency 
range.  Although  this  range  was  extended  down  to  60  kHz  for  the  laboratory  experiments,  it 
is  still  not  low  enough  to  observe  significant  dispersion  for  glass  beads  saturated  with  water. 
To  overcome  this  limitation,  a  second  glass  bead  sediment  was  saturated  with  Silicone  oil 
which  has  a  viscosity  of  approximately  0.1  Pa*s.  Biot  theory  predicts  that  by  increasing  the 
viscosity,  the  region  of  significant  dispersion  should  shift  upward  in  frequency  and  hence 
should  be  observable  with  the  measurement  system.  Details  of  the  sediments  as  well  as 
the  measurement  system  are  discussed  in  Section  II.  This  experiment  is  similar  to  the 
measurements  made  by  Seifert  et.al.13  In  those  measurements,  several  different  fluids  were 
used  to  vary  the  viscosity  of  the  pore  fluid  from  0.001  to  1  Pa-s.  However,  they  observed  an 
f2  frequency  dependence  in  the  attenuation  whereas  the  measurements  presented  here  show 
an  Z1  dependence  over  a  similar  range  of  frequencies.  For  the  measurements  in  a  silicone 
oil  saturated  sediment,  significant  dispersion  is  observed  in  the  current  experiments.  These 
results  are  similar  to  those  of  SAX99  in  that  the  attenuation  deviates  from  Biot  theory 
over  the  same  frequency  range  while  the  sound  speed  is  well  described  by  Biot  theory.  The 
measurement  results  and  model  comparisons  are  shown  in  Section  III.  These  results  do 
not  support  the  scattering  hypothesis  for  the  SAX99  results  and  while  neither  Biot  nor 
Buckingham’s  theory  individually  captures  all  aspects  of  the  data,  a  model  which  combines 
the  important  aspects  of  both  is  able  to  describe  all  three  data  sets.  This  model  is  developed 
in  Section  IV  while  the  implications  of  the  data-model  comparisons  using  this  hybrid  model 
are  discussed  in  Section  V. 

II.  EXPERIMENTAL  PROCEDURE 

The  system  used  to  make  the  measurements  described  here  was  deployed  during  SAX99 
by  APL-UW.10’12  This  diver-deployed  system  was  known  as  the  ‘attenuation  array’  and 
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was  able  to  make  sound  speed  measurements  in  the  100-200  kHz  frequency  range  and 
attenuation  measurements  in  the  100-260  kHz  frequency  range.  The  system  consists  of  four 
ITC  6148  transducers  (two  receivers  and  two  transmitters)  which  allowed  measurements  to 
be  made  over  four  separate  paths.  During  the  experiment  the  path  lengths  were  from  22  to 
44  cm  and  the  transducers  were  inserted  into  the  sediment  to  a  depth  of  10  cm.  Because 
this  system  was  designed  to  be  diver-deployed,  it  was  relatively  compact  and  thus  easy  to 
use  for  laboratory  based  measurements. 

The  sediments  in  the  laboratory  measurements  were  composed  of  well-sorted  glass  beads 
which  had  diameters  in  the  range  of  0.43  to  0.30  mm  and  a  minimum  roundness  of  90%. 
Two  fluids  were  used  to  saturate  the  glass  beads:  water  and  a  silicone  oil  with  a  viscosity 
of  0.098  Pa-s.  Biot  theory  predicts  that  by  increasing  the  viscosity  of  the  pore  fluid,  the 
dispersion  in  the  sound  speed  should  shift  to  higher  frequencies7.  By  shifting  the  dispersion 
up  in  frequency,  it  should  be  observable  in  the  frequency  range  of  the  attenuation  array.  It 
should  also  provide  a  second  test  of  Biot  theory  in  that  the  attenuation  should  scale  with 
the  viscosity  of  the  pore  fluid.  In  the  measurements  of  Seifert  et.  al.  this  scaling  was  not 
observed  leading  them  to  the  reject  Biot’s  porous  medium  description  of  the  unconsolidated 
sand  sediment13. 

In  order  to  insure  that  bubbles  would  not  be  present  in  the  sediments,  a  vacuum  was 
applied  to  the  fluid  before  the  sediment  was  added  and  allowed  to  degas  over  several  days. 
In  order  to  withstand  the  vacuum,  two  stainless  steel  oil  drums  were  chosen  as  tanks 
for  the  sediments.  These  drums  had  an  inner  diameter  of  54.6  cm  and  a  depth  of  83.5 
cm.  The  glass  beads  were  added  to  the  degassed  fluids  in  10  cm  layers  and  the  glass 
bead/fluid  combination  was  degassed  and  then  stirred  after  each  layer  was  added.  Over 
the  course  of  the  experiments,  the  sediments  were  stirred  and  a  vacuum  was  applied  to  the 
tanks  periodically  to  remove  any  bubbles  that  may  have  formed.  The  glass  bead/water 
sediment  had  a  depth  of  37  cm  with  a  water  layer  above  it  that  was  36  cm  thick.  The  glass 
bead/Silicone  oil  sediment  had  a  depth  of  23  cm  with  a  22  cm  thick  layer  of  Silicone  oil 
above  it. 

The  properties  of  each  sediment,  including  that  of  SAX99,  are  given  in  Table  1.  The 
bulk  modulus  and  density  of  the  glass  beads  were  provided  by  the  manufacturer,  Potters 
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Industries  Inc.,  and  are  consistent  with  reported  values  of  glass  found  in  the  literature.  The 
bulk  modulus  of  the  water  was  determined  from  the  average  sound  speed  measurement  in 
the  water.  The  density  and  viscosity  of  the  Silicone  oil  was  provided  by  the  manufacturer, 
Silchem,  and  the  bulk  modulus  of  the  oil  was  determined  from  its  measured  sound  speed. 

The  remaining  properties  important  for  the  application  of  Biot  theory  and  which  are 
often  difficult  to  measure  in  sand  sediments  are  the  porosity,  /?,  permeability,  k,  and  the 
tortuosity,  a.  The  permeability  of  the  glass  bead/water  sediment  was  determined  by  taking 
samples  of  the  sediment  with  a  known  volume  and  measuring  the  density  before  and  after 
removing  the  water  by  drying  the  sample  in  an  oven.  The  value  of  0.384  is  reasonably  close 
to  the  random  close  packing  limit  for  identical  spheres.  The  same  procedure  could  not  be 
applied  to  the  Silicone  oil  sediment  because  the  vapor  pressure  of  the  oil  is  much  higher 
than  that  of  water.  The  random  close  packed  limit  of  0.36  was  assumed  here  instead.  It 
may  seem  reasonable  to  assume  that  the  porosity  should  be  the  same  in  both  cases,  but 
this  assumption  was  rejected  because  the  macroscopic  properties  of  the  sediments  were 
found  to  be  much  different.  The  Silicone  oil  sediment  was  found  to  heal  much  faster  than 
the  water  sediment  when  stirred  or  when  transducers  were  inserted  and  removed  from  the 
sediment.  Also,  during  degassing,  bubbles  were  able  to  move  out  from  the  interior  of  the 
silicone  oil  sediment  much  faster  and  easier  than  from  the  water  saturated  sediment.  It 
is  hypothesized  that  this  is  due  to  the  higher  viscosity  of  the  silicone  oil  which  allows  the 
beads  to  slide  relative  to  one  another  with  less  friction  due  to  hydrodynamic  lubrication14. 
It  seems  more  likely  that  the  beads  are  able  to  reach  a  denser  ground  state  in  Silicone  oil 
than  they  would  be  able  to  in  water. 

The  hydraulic  permeability  was  calculated  using  the  Kozeny-Carmen  relation15, 


180  (1  —  /?)2’  ^ 
where  d  is  the  average  diameter  of  the  glass  beads.  This  relationship  was  tested  for  the 
water /glass  bead  sediment  using  a  constant  head  permeameter  and  the  value  calculated 
using  the  Kozeny-Carmen  relation  was  found  to  fall  within  the  error  of  the  measurement.  A 
measurement  of  the  tortuosity  is  notoriously  difficult  to  make  for  unconsolidated  sediment 
and  as  a  result  an  expression  derived  by  Berryman  for  a  collection  of  spheres  was  used16. 
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The  expression, 


2 

a“  1  +  /?’ 

relates  the  tortuosity  to  the  porosity  and  gives  a  lower  value  than  the  expression  that  was 
used  previously  by  Berryman  to  describe  sintered  glass  bead  packs17.  As  Berryman  shows 
in  Table  1  of  Ref.  16,  the  expression  given  in  Eq.  2  fits  the  analytical  values  for  ordered 
unconsolidated  glass  bead  packs  whereas  the  previous  expression  is  a  better  fit  to  sintered 
glass  bead  packs.  As  a  result,  Eq.  2  is  used  for  the  tortuosities  of  the  glass  bead  sediments. 

Because  the  diameter  of  the  tanks  were  smaller  than  the  original  attenuation  array,  the 
array  was  reconfigured  for  use  in  the  tanks.  The  new  arrangement  maintained  the  original 
electronics  and  transducer  mounts,  while  repositioning  the  transducers  such  that  they  sat 
at  the  corners  of  a  rectangle.  For  each  transmitter,  one  receiver  was  mounted  close  to 
provide  a  short  path  length,  while  the  second  transducer  was  mounted  further  away  to 
provide  a  longer  path  length.  The  sound  speed  was  determined  using  the  difference  in 
phases  between  a  pulse  received  at  two  different  transducers18.  As  a  result,  the  difference 
in  path  lengths,  A P,  is  the  important  distance  for  this  measurement.  For  the  water/glass 
bead  sediment  A  Pi  =  25.99  cm  and  AP2  ==  25.95  cm  while  for  the  silicone  oil/glass  bead 
sediment  A  Pi  =  14.39  cm  and  AP2  —  14.97  cm.  The  reduction  in  path  length  for  the  silicone 
oil  was  necessary  in  order  to  overcome  the  signal  loss  due  to  the  increased  attenuation.  The 
attenuation  was  determined  from  the  spectral  amplitudes  using  the  method  described  in 
Ref.  18.  This  technique  also  depends  on  the  difference  in  the  two  path  lengths.  The  sound 
speed  and  attenuation  measurements  for  both  sediments  were  performed  using  the  original 
attenuation  array  over  a  frequency  range  of  100  kHz  to  260  kHz.  However,  due  to  the  large 
attenuation,  measurements  above  200  kHz  could  not  be  made  in  the  silicone  oil  saturated 
sediment. 

For  each  measurement,  the  transducers  were  first  placed  in  the  overlying  fluid  layer 
and  the  sound  speed  and  spectral  amplitude  of  the  pulse  in  the  fluid  were  measured.  The 
transducers  were  then  inserted  into  the  sediment.  The  depth  of  the  tranducers  in  the 
water /glass  bead  sediment  was  10  cm  while  in  the  silicone  oil/glass  bead  sediment  it  was 
11.4  cm.  The  tank  was  then  shaken  by  hand  for  approximately  two  minutes  then  allowed 
to  sit  for  4-6  hours  to  allow  the  sediment  to  settle.  The  sound  speed  and  spectral  amplitude 
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were  then  measured  and  after  the  transducers  were  removed  from  the  sediment,  the  sound 
speed  and  spectral  amplitudes  in  the  fluid  were  again  measured.  Once  the  transducers  were 
removed  from  the  tank,  the  sediment  was  again  shaken  and  allowed  to  settle  overnight. 

In  order  to  extend  the  lower  limit  of  the  frequency  range  of  the  array,  a  Reson  TC4040 
was  introduced  which  was  able  to  operate  in  the  60  kHz  to  100  kHz  range.  The  difference 
in  path  length  for  the  new  array  was  AP\  =  7.82  cm  in  water  and  AP\  =  4.53  cm  in  silicone 
oil.  Unfortunately,  there  was  an  unresolved  problem  with  the  phase  of  the  output  pulse  of 
the  new  array  which  produced  an  offset  in  the  measured  sound  speed.  Consequently  these 
phase  speed  measurements  were  discarded.  The  attenuation  measurements  were  unaffected 
by  this  problem  and  are  presented  here. 

III.  RESULTS  AND  MODEL-DATA  COMPARISONS 

The  sound  speed  ratio  and  attenuation  measured  in  each  sediment  are  presented  in 
Fig.  2.  Each  data  point  is  the  mean  of  measurements  of  approximately  20  realizations 
of  the  each  sediment.  The  error  bars  on  each  data  point  reflect  the  uncertainties  in  the 
path  lengths  as  well  the  distribution  of  the  measured  values.  As  predicted  by  Biot  theory, 
the  sound  speed  in  the  silicone  oil/glass  bead  sediment  exhibits  greater  dispersion  then  the 
water/glass  bead  sediment.  However,  contrary  to  Biot  theory,  both  media  show  a  frequency 
dependence  for  the  attenuation  that  is  greater  than  / 1/2.  A  linear  fit  to  the  attenuation 
shows  that  the  attenuation  follows  an  fim  dependence  for  the  water/glass  bead  sediment 
and  an  /°  "  dependence  for  the  silicone  oil/glass  bead  sediment. 

In  order  to  test  Biot  theory  against  this  data,  it  is  necessary  to  know  certain  material 
parameters  in  order  to  determine  the  sound  speed  and  attenuation.  A  majority  of  these 
parameters  are  discussed  in  the  previous  section  and  are  given  in  Table  I,  but  there  are 
two  notable  exceptions.  For  both  glass  bead  sediments,  the  frame  bulk  modulus,  A&,  and 
the  shear  modulus,  //,  are  unknown.  Williams  has  shown  that  in  most  sand  sediments, 
because  the  bulk  and  shear  frame  moduli  are  much  smaller  than  the  other  moduli,  they 
can  be  neglected  with  only  minor  errors  (0.5%)  in  the  high  frequency  sound  speed  and 
attenuation.  Taking  the  Kb  —  fi  =  0  limit  of  Biot  theory  yields  the  effective  density  fluid 
model  (EDFM)  which  gives  the  complex  wavenumber  for  the  compressional  wave.  The 
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EDFM  along  with  the  fluid,  solid,  and  sediment  properties  in  Table  I  yield  the  solid  lines 
in  Fig.  2.  In  both  sediments,  the  EDFM  captures  the  measured  sound  speed  values  but 
does  not  capture  the  high  frequency  behavior  of  the  attenuation.  Although  the  EDFM 
does  not  capture  the  details  of  the  attenuation,  it  does  predict  an  increase  in  attenuation 
between  the  water  and  silicone  oil  sediments  that  is  on  the  same  order  as  the  increase  in  the 
measured  attenuation.  Consequently,  both  the  frequency  scaling  and  frequency  dependence 
contradict  what  was  measured  by  Seifert  et.  al.13 

The  compressional  wavenumber  in  Buckingham’s  model  can  be  expressed  as4 

2  uPp 

H  +  It  (—iu>T)n  ’  1  j 

where  u  is  the  radial  frequency,  p  =  (1  —  f$)ps  +  j3pf  is  the  average  density  of  the  sediment, 
H  =  (fi/Kf  +  (1  —  j3)/Ks)~l  is  the  bulk  modulus  of  the  sediment  evaluated  from  Wood’s 
equation,  7t  =  7P  +  4/37 s  the  rigidity  coefficient,  and  m  is  the  material  exponent.  Note 
that  the  rigidity  coefficient  is  a  combination  of  Buckingham’s  compressional  and  shear 
rigidity  coefficients.  Since  there  is  not  a  shear  wave  measurement  to  use  to  determine  the 
values  of  m  and  7S,  the  rigidity  coefficients  can  be  combined  into  a  single  constant  which 
must  be  fit  to  the  data  along  with  the  material  exponent.  Note  that  the  assumption  of 
m  —  n  has  also  been  made. 

The  sound  speed  and  attenuation  given  by  Buckingham’s  model  are  shown  as  the  dashed 
lines  in  Fig.  2  and  the  values  of  and  m  are  given  in  Table  I.  For  the  water /glass  bead 
sediment,  the  values  of  these  parameters  are  comparable  to  those  determined  for  the  SAX99 
data  in  that  m  and  7 t  are  the  same  order  of  magnitude.  For  the  silicone  oil  sediment,  m  is 
an  order  of  magnitude  larger  while  7 t  is  two  orders  of  magnitude  smaller.  It  may  be  that 
the  large  viscosity  of  the  silicone  oil  brings  into  doubt  the  approximation  made  in  Eq.  58 
and  59  in  Ref.  4  which  stems  from  the  assumption  that  viscous  effects  at  the  grain  contacts 
are  negligible.  If  these  approximations  are  not  made,  the  number  of  unknowns  increases 
to  three  (this  again  assumes  that  the  shear  and  compressional  values  are  the  same),  the 
determination  of  which  is  beyond  the  scope  of  this  paper. 
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the  components  of  the  strain  tensor  and  how  they  relate  to  the  stresses  given  by  Eq.  (4) 
and  (5).  Consider  first  the  shear  strain  in  the  frame  matrix, 


1 


l  1  (dui  duj  . 

-»  27"  =  2fe  +  Slj'  (6) 

This  strain  will  produce  a  shear  stress  within  the  matrix  due  to  the  translational  sliding  of 
the  grains  against  on  another.  From  Eq.  (5),  the  shear  stress  is 


(7) 


where  Vi  =  it*.  Assuming  that  the  disturbance  in  the  matrix  has  a  harmonic  time  depen¬ 
dence,  uie~twt ,  then  Eq.  (7)  becomes 


cr,j  =  2fib  (u>)  €ij 


(8) 


where 


(xh  (u)  =  -iiur)sH*  (ioS)  (9) 

and  H*  ( iu ;)  is  the  complex  conjugate  of  the  temporal  fourier  transform  of  hs  ( t ). 

For  the  second  elastic  modulus  of  the  granular  matrix,  consider  the  response  of  a  unit 
cube  of  the  matrix  to  a  uniform  pressure.  This  pressure  gives  rise  to  a  change  in  the  volume 
of  the  cube  equal  to  e  =  Vw  and  will  also  produce  radial  sliding  at  the  grain  contacts  due 
to  the  compression  of  the  grains.  Using  Eq.  (4),  the  stress-strain  relation  in  this  situation 
becomes 


cr  =  APV  [hp  (t)  ®  v\ .  (10) 

Once  again  assuming  a  harmonic  compression,  the  stress-strain  relation  becomes 

a^Kb(u)e,  (11) 

where 

Kb  (u)  =  —iujXpH*  (iu)  (12) 
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and  is  the  complex  conjugate  of  the  temporal  fourier  transform  of  hp(t).  The 

dynamic  elastic  moduli  of  the  matrix  are  therefore  the  grain-to-grain  shear  modulus,  (a;), 
and  the  grain-to-grain  bulk  modulus,  Kb  (a;). 

The  temporal  Fourier  transforms  of  the  MIRFS  in  Eq.  (9)  and  (12)  are  given  approxi¬ 
mately  by4 


and 


h;m 


r(l  -m) 


(13) 


F  (1  —  m) 


(14) 


where  Y  is  the  complement  of  the  incomplete  gamma  function,  m,  n  are  the  material  ex¬ 
ponents,  and  tSiP  are  the  stress  relaxation  time  constants.  The  approximations  are  valid 
assuming  that  ojts#  1.  With  these  expressions,  the  elastic  moduli  can  be  written  as 


Vb(oj)  =  7 5  {—iojT)m 

and 

Kh(u)  =  lp(-iojT)n 

where  the  shear  and  compressional  rigidity  coefficients  are 

(I-)©’  <-> 

and 

The  time  constant  T  has  been  introduced  to  keep  terms  raised  to  the  fractional  powers  n 
and  m  dimensionless.  The  expressions  are  independent  of  T  and  hence  it  is  conveniently 
taken  to  be  T  =  Is.  These  expressions  for  the  frame  moduli  can  be  used  in  Biot’s  equations. 
A  detailed  discussion  of  Biot’s  equations  can  be  found  in  Ref.  9  while  a  concise  presentation 
is  found  in  Eq.  (l)-(8)  of  Ref.  12.  The  use  of  Eq.  (15)  and  (16)  in  Biot  Theory  will  be 
referred  to  as  the  “Biot- Buckingham  model”  for  the  remainder  of  this  paper. 


(15) 

(16) 
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The  expressions  for  the  moduli  given  by  Eq.  (15)  and  (16)  are  not  new  and  were  first 
suggested  by  Biot  in  the  development  of  a  phenomenological  approach  to  solid  dissipation 
in  porous  medium6,7  and  later  explored  in  detail,  outside  the  context  of  Biot  theory,  by 
Kjartansson20.  Kjartansson  was  motivated  by  early  measurements  on  absorption  of  sound 
in  rocks  which  found  that  the  energy  lost  per  cycle  or  wavelength,  Q,  was  independent  of 
the  frequency.  To  account  for  this,  he  proposed  a  model  for  the  elastic  moduli  which  has 
essentially  the  same  form  as  Eq.(15)  and  (16).  For  the  shear  modulus,  and  likewise  for 
the  bulk  modulus,  the  phase  angle  between  the  stress  and  the  strain  can  be  determined  by 
rewriting  the  modulus  as 


Hb  (oj)  =  7*  |wT|m  exp  (— *-7j-sgn  (w))  .  (19) 

The  Q  is  directly  related  to  this  phase  angle  and  from  this  expression  for  the  shear  modulus, 


1  (m  7r> 

—  =  tan  i  - 


Q=,“V2 


-)  •  (20) 

Because  of  this  property,  the  model  is  often  referred  to  as  the  constant  Q  (CQ)  model. 

Kjartansson  showed  that  this  model  could  be  expressed  as  an  infinite  set  of  viscoelastic 
Zener  elements20.  As  a  consequence  of  this  construction,  the  CQ  model  is  causal.  Buck¬ 
ingham  has  also  shown  this  directly  in  the  development  of  his  model2 .  The  introduction  of 
CQ  frame  moduli  into  Biot  theory  was  investigated  by  Turget21  who  numerically  integrated 
the  Kramers-Kronig  relations  and  also  found  that  this  viscoelastic  Biot  theory  was  causal. 

To  examine  the  properties  of  this  hybrid  model,  it  is  easiest  to  look  first  at  the  behavior 
of  the  shear  wave.  In  the  Biot-Buckingham  model,  the  shear  wave  number  becomes 


:  U) 


(. s ) 

Peff 

M") 


(21) 


where 


Peff 


Pf 


'a  (1  —  /?)  ps  +  f3  (a  —  1)  + 


ippyFM ' 

PfKLJ 


OL  + 


ipijFfa) 


(22) 


and  (it  is  the  viscoelastic  shear  modulus  given  by  Eq.  (15).  For  a  detailed  derivation  of 
this  expression  see  p.  19  of  Ref.  9.  All  of  the  constants  in  p^jj  can  be  measured  directly. 
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Following  Buckingham,  the  material  exponent,  ra,  and  the  rigidity  coefficient,  7,,  can  be 
evaluated  given  a  measurement  of  the  shear  speed  and  attenuation  at  a  single  frequency. 
Writing  the  shear  modulus  as 


7S  (— iioT)m  =  7S  (c uT)m  (cos  (to7t/2)  —  i  sin  (m7r/2)) ,  (23) 

the  material  exponent  can  then  be  found  from 


2  lm  (Peff/k2s) 


and  the  rigidity  coefficient  is  then 


(24) 


V*Peff  ^ 

7s  (—iu)T)m  k2s  ’  { ’ 

where  ks  =  u>/cs  —  cs  is  the  shear  wave  speed,  and  as  is  the  attenuation.  During 
SAX99,  the  shear  wave  speed  and  attenuation  were  measured11  at  1  kHz  and  were  cs  =  120 
m/s  with  a  range  of  97  to  147  m/s  and  as  =  30  dB/m  with  a  range  of  21  to  40  dB/m. 
Using  these  values,  the  material  exponent,  m,  has  a  best  fit  of  0.0506  to  the  mean  sound 
speed  and  attenuation.  The  rigidity  coefficient  has  a  best  fit  of  1.82  xlO7  Pa  to  the  mean 
values  of  the  sound  speed  and  attenuation. 

The  shear  speed  and  attenuation  are  shown  as  the  solid  lines  in  Fig.  3.  Also  shown  in  this 
figure  are  the  predictions  of  the  Buckingham  model  and  the  predictions  of  Biot  theory  with 
a  complex  shear  modulus  with  no  frequency  dependence.  For  this  complex  shear  modulus, 
the  best  fit  to  the  mean  values  of  the  sound  speed  and  attenuation  is  2.83  x  107  —  £2.25  x  106 
Pa. 


In  Ref.  2,  Buckingham  attempts  to  address  the  question  of  whether  it  is  necessary  to 
account  for  viscous  losses  due  to  global  of  local  motion  of  the  pore  fluid  when  modeling 
sound  propagation  in  marine  sediments.  In  order  to  do  this,  he  introduces  viscous  loss 
terms  into  the  equation  of  motion  for  the  saturated  sediment.  The  solution  to  the  resulting 
wave  equation  indicates  a  frequency  dependence  for  both  the  sound  speed  and  attenuation 
from  which  he  is  able  to  draw  two  conclusions.  In  the  combined  model,  the  intergranular 
friction  appears  to  dominate  at  lower  frequencies  while  the  pore  fluid  viscosity  dominates  at 
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high  frequencies.  Also,  the  sound  speed  varies  with  frequency  as  f1/2  in  the  high  frequency 
regime  where  the  pore  fluid  viscosity  dominates. 

In  the  Biot-Buckingham  theory,  the  effects  of  both  loss  mechanisms  are  not  as  clear  cut  as 
Buckingham’s  test  model.  In  the  low  frequency  limit,  the  sound  speed  varies  with  frequency 
as  /m/2,  while  the  attenuation  varies  as  /2-m/2  for  m  «  1  and  then  transitions  to  / 1~rn/2 
as  m  increases.  For  the  properties  of  the  SAX99  sediment,  which  are  typical  of  a  sandy 
sediment,  this  transition  occurs  for  extremely  small  m  and  by  m  «  0.02  the  intergranular 
friction  has  become  the  dominant  loss  mechanism,  as  predicted  by  Buckingham’s  test  model. 
In  the  high  frequency  limit,  the  shear  speed  also  varies  with  frequency  as  /m/2.  This 
contradicts  Buckingham’s  test  model  where  the  sound  speed  varied  as  Z1/2,  indicating  that 
even  at  high  frequencies  the  intergranular  friction  remains  dominate.  For  high  frequencies, 
the  attenuation  also  exhibits  behavior  similar  to  what  was  observed  at  low  frequencies,  for 
small  m,  the  attenuation  varies  with  frequency  as  J1/2-™/2,  but  once  again  as  soon  asm  >  0 
the  frequency  dependence  begins  a  transition  to  /1_m/2  which  is  complete  by  m  «  0.02. 

The  asymptotics  of  Eq.  (21)  indicate  that  in  both  high  and  low  frequency  limits,  the 
intergranular  friction  is  the  dominant  loss  mechanism,  while  from  Fig.  3  both  mechanisms 
play  important  roles  in  the  mid-frequency  range.  This  is  similar  to  the  results  of  Biot- 
Stoll  theory.  Brunson  has  shown  that  Biot  theory  with  complex  shear  moduli  is  able  to 
account  for  details  of  the  measured  attenuation,  specifically  the  “knee”  in  Fig.  6  of  Ref.  22 
that  a  linear  fit  would  miss.  (See  also  Fig.  1  of  Ref.  3  for  further  deviations  from  the 
/ 1  dependence  for  shear  attenuation.)  From  the  similarities  in  the  behavior  of  the  Biot- 
Buckingham  and  Biot-Stoll  theories  in  Fig.  3  it  seems  that  Biot-Buckingham  theory  would 
also  be  able  to  capture  these  details  in  the  attenuation.  A  true  test,  as  seen  in  Fig.  3(a) 
would  be  a  measurement  of  the  shear  speed  across  a  large  frequency  range,  however  this 
measurement  may  be  prohibitively  difficult  in  sandy  sediments. 

Using  the  values  for  the  shear  rigidity  coefficient  and  the  corresponding  material  expo¬ 
nent,  the  remaining  coefficient  and  exponent  can  be  determined  for  the  compressional  wave 
through  a  best  fit  to  the  data.  This  approach  differs  from  Buckingham’s  development  of  his 
model  where  he  assumes  that  the  material  exponents  should  be  equal  to  one  another4.  In 
his  model,  the  grain-to-grain  friction  is  due  to  a  number  of  random  micro-asperity  sliding 
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events,  each  of  one  which  is  assumed  to  be  identical  regardless  of  whether  the  grain-to- 
grain  motion  is  radial  or  compressional.  Since  the  material  exponents  characterize  these 
individual  events,  it  seems  reasonable  that  they  should  have  the  same  values.  Buckingham 
presents  limited  evidence  in  the  form  of  data-model  comparisons  that  support  this  con¬ 
clusion  for  his  model.  The  data-model  comparisons  presented  here  support  the  conclusion 
that  n  and  m  are  not  equal,  hence  they  are  treated  separately. 

The  compressional  sound  speed  and  attenuation  predicted  by  the  Biot-Buckingham 
model  for  the  SAX99  sediment  are  presented  in  Fig.  4(a)  and  (d).  Both  n  and  were 
determined  from  a  best  fit  to  the  data  and  were  found  to  be  n  =  1.088  and  7 p  =  41.96  Pa. 
This  hybrid  model  succeeds  in  combining  the  strengths  of  both  the  Biot  and  Buckingham 
models  and  follows  the  data  better  than  either  individual  model.  The  model  also  contra¬ 
dicts  Buckingham’s  conclusions  about  the  importance  of  losses  due  to  viscous  flow  for  the 
compressional  wave.  Whereas  he  concludes  that  viscous  loss  is  the  dominant  mechanism 
at  high  frequencies,  the  Biot-Buckingham  model  predicts  that  viscous  loss  dominates  in 
the  mid-frequencies  while  granular  friction  is  the  dominant  mechanism  at  high  and  low 
frequencies.  This  is  similar  to  the  asymptotics  of  the  shear  wave  where  granular  friction 
also  dominates  at  both  the  high  and  low  frequencies.  An  analysis  of  the  asymptotics  of 
the  compressional  wave  are  not  as  straightforward  as  for  the  shear  wave,  but  there  are 
several  observations  which  can  be  made  immediately.  Since  the  frequency  dependence  of 
the  frame  bulk  modulus  is  nearly  linear,  for  low  frequencies  granular  friction  will  dominate 
the  attenuation  since  the  attenuation  due  to  viscous  flow  goes  as  /2.  In  the  high  frequency 
limit,  Kb  (a;)  becomes  very  large  and  hence  dominates  the  wave  propagation. 

In  order  to  compare  the  Biot-Buckingham  model  to  the  data  from  the  glass  bead  sed¬ 
iments,  four  variables  need  to  be  determined:  n,  m,  7P,  and  7*.  Unlike  the  SAX99  data 
there  are  no  shear  wave  measurements  to  use  to  determine  the  shear  rigidity  coefficient 
and  exponent.  The  results  for  the  SAX99  data,  however,  suggest  that  the  shear  modulus 
may  be  negligible  compared  to  the  other  moduli.  Setting  fib  =  0  for  the  Biot-Buckingham 
model  yields  both  sound  speed  and  attenuation  predictions  that  are  nearly  indistinguish¬ 
able  from  the  EDFM  predictions.  Consequently,  for  the  glass  bead  sediments,  the  shear 
rigidity  coefficient  is  taken  as  7S  =  0.  The  predictions  of  the  Biot-Buckingham  model  are 
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shown  in  Fig.  4(b)  and  (e)  for  the  water  and  glass  bead  sediment  and  in  Fig.  4(c)  and  (f) 
for  the  silicone  oil  and  glass  bead  sediments. 

For  water  and  glass  beads,  the  grain- to-grain  shearing  parameters  were  determined  to 
be  n  =  0.915  and  7P  =  182.94  Pa.  The  predictions  match  the  data  well  in  this  case, 
but  it  is  difficult  to  say  whether  or  not  is  does  any  better  than  the  Buckingham  model 
over  this  limited  frequency  range.  For  the  silicone  oil  and  glass  beads,  the  grain-to-grain 
shearing  parameters  were  determined  to  be  n  =  1.059  and  7P  =  57.89  Pa.  The  Biot- 
Buckingham  model  captures  the  data  very  well  in  this  case  and  with  values  for  the  grain- 
to-grain  parameters  that  are  consistent  with  the  other  sediments.  This  differs  from  the 
Buckingham  model  for  which  the  material  exponent  increases  by  a  factor  of  10  relative  the 
water/glass  bead  sediment.  This  indicates  that  for  the  Buckingham  model  the  increase 
in  attenuation  and  dispersion  must  be  due  to  a  decrease  in  the  strain-hardening  at  the 
micro-asperities  even  though  the  viscosity  of  the  pore  fluid  was  increased  by  a  factor  of 
100.  In  the  Biot-Buckingham  model,  the  increase  in  attenuation  is  due  to  the  increase  in 
the  viscosity  of  the  pore  fluid,  while  the  attenuation  due  to  the  grain-to-grain  friction  is 
nearly  identical  to  that  observed  in  the  previous  glass  bead  sediment  and  in  the  SAX99 
sediment. 

Not  only  are  the  grain-to-grain  shearing  parameters  similar  in  all  three  sediments,  but 
the  material  exponent  is  very  nearly  equal  to  one,  as  can  be  seen  in  Table  II.  This  suggests 
that  the  Biot-Buckingham  model  with  n  =  1  may  yield  predictions  that  compare  favorably 
to  the  data.  The  results  of  the  Biot-Buckingham  model  with  n  =  1  and  7p  taken  as  the  free 
parameter,  are  shown  as  the  dashed  lines  in  Fig.  4.  For  the  attenuation,  the  two  curves  are 
practically  indistinguishable  while  for  the  sound  speeds  the  deviation  is  most  pronounced 
at  high  frequencies  where  the  grain-to-grain  losses  are  greatest.  Even  with  this  deviation, 
the  model  still  compares  favorably  to  the  data. 

V.  DISCUSSION 

The  success  of  the  Biot-Buckingham  model  with  n  =  1  has  implications  for  the  theo¬ 
retical  framework  of  the  grain-to-grain  shearing  mechanism  that  Buckingham  develops  in 
Ref.  4.  The  material  exponent  is  defined  as  n  =  EP/6P,  where  Ep  is  the  elastic  response 
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of  the  micro-asperity  contact  and  6p  reflects  the  nonlinear  response  of  the  viscous  damping 
at  the  contact.  Thus  n  =  1  means  that  elastic  and  viscous  response  is  equal.  This  un¬ 
likely  equality  brings  Buckingham’s  theoretical  framework  into  question  and  suggests  that 
another  explanation  for  the  form  of  Eq.  16  may  be  needed. 

From  a  purely  phenomenological  standpoint,  the  complex  moduli  can  be  constructed 
from  simple  elastic  and  viscous  elements.  Consider  the  form  of  the  bulk  modulus  in  the 
Biot-Buckingham  model  when  n  =  1, 

Kb  (c o)  =  (26) 

The  modulus  is  purely  imaginary  and  corresponds  to  a  viscoelastic  model  composed  entirely 
of  dashpots.  In  Biot’s  phenomenological  viscoelastic  model,  he  suggests  a  viscoelastic 
modulus  of  the  form, 

Kb(co)  =  a  (—iu>)s  +  Kb,  (27) 

which,  for  s  =  1,  corresponds  to  the  modulus  in  Eq.  26  in  series  with  a  real  elastic  modulus. 
For  small  values  of  this  elastic  term,  the  results  of  using  Eq.  27  will  not  differ  significantly 
from  the  predictions  shown  in  Fig.  4. 

Biot  suggests  several  other  expressions  for  the  viscoelastic  behavior  of  the  moduli  and 
suggests  some  possible  mechanisms  which  might  lead  to  viscoelasticity  in  porous  media. 
Among  these  mechanisms,  he  suggests  that  dissipation  may  be  present  at  the  granular 
contacts  due  the  viscous  flow  of  fluid  being  squeezed  in  and  out  of  the  contact  region. 
Biot  constructs  a  simple  spring  and  dashpot  model  of  this  squeeze  flow  but  he  is  quick 
to  recognize  that  this  simple  model  is  deficient.  The  viscous  response  of  the  fluid  at  the 
scales  of  the  intergranular  gaps  differs  considerably  from  the  response  of  the  bulk  fluid. 
Subsequent  work  in  tribology  supports  this  conclusion  and  Buckingham’s  model  is  the  first 
to  try  to  capture  this  complex  fluid-solid  interaction  at  the  grain  contacts. 

Although  it  is  not  possible  to  draw  any  conclusions  from  the  theoretical  implications  of 
the  data-model  comparisons,  the  fits  of  the  Biot-Buckingham  model  do  provide  some  insight 
into  the  nature  of  the  attenuation  mechanisms.  The  increase  in  the  attenuation  when  the 
pore  fluid  is  changed  from  water  to  silicone  oil  is  captured  by  Biot  theory  indicating  that 
viscous  flow  is  an  important  contributor  to  the  acoustic  loss.  For  a  viscosity  increase  of 
two  orders  of  magnitude,  the  rigidity  coefficient  doubles  as  seen  in  Table  II  for  n  =  1. 
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The  deviation  of  the  rigidity  coefficient  for  the  SAX99  sediment  from  this  simple  viscosity 
dependence  suggests  that  the  attenuation  mechanism  may  depend  in  a  complex  way  on  the 
other  properties  of  the  sand  grain  such  as  the  surface  roughness  of  the  grain. 

The  discussion  here  has  focused  mainly  on  intergranular  dissipation  as  the  possible 
explanation  for  the  deviation  of  the  attenuation  from  the  predictions  of  Biot  theory.  The 
initial  hypothesis  was  that  scattering  was  responsible  for  the  losses  observed  in  the  SAX99 
data.  Although  the  measurements  in  glass  beads  suggest  that  the  excess  attenuation  is  not 
due  to  scattering,  there  is  still  evidence  of  scattering  in  the  SAX99  data  that  is  remains 
to  be  addressed.  The  transmitted  waveforms  in  the  sediment  show  significant  distortion 
and  degradation  as  they  propagate  and  potential  scatterers  are  known  to  be  present  in  the 
sediment12.  This  kind  of  signal  degradation  was  not  observed  in  the  glass  bead  sediments 
suggesting  that  scattering  may  still  be  a  contributor  to  the  acoustic  attenuation  in  real 
sediments.  The  measurements  presented  here  do  suggest,  though,  that  it  may  not  be 
the  sole  or  dominant  attenuation  mechanism  at  high  frequencies  and  that  a  mechanism 
more  fundamental  to  the  nature  of  unconsolidated  porous  media  needs  to  be  taken  into 
account.  Whether  these  losses  are  related  to  grain-to-grain  interactions  or  to  some  other 
mechanism,  such  as  acoustic  interactions  with  “force  chains,”  remains  an  open  question  for 
further  research  and  speculation. 
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TABLE  I:  Parameters  used  in  calculating  the  effective  density  fluid  model  and  Buckingham  model 
results.  The  Biot-Buckingham  model  results  were  also  calculated  using  these  results  except  for 
the  material  exponent  and  rigidity  coefficient. 


Parameter 

Value 

SAX99 

Water/Beads 

SilOO/Beads 

Bulk  modulus  of  grains  or  glass  beads  ( Kg ) 

3.2  x  1010  Pa 

7.0  x  1010  Pa 

7.0  x  1010  Pa 

Mass  density  of  grains  or  glass  beads  ( pg ) 

2690  kg/m3 

2500  kg/m3 

2500  kg/m3 

Bulk  modulus  of  pore  fluid  ( Kf ) 

2.395  x  109  Pa 

2.23  x  109  Pa 

9.68  x  108  Pa 

Mass  density  of  pore  fluid  (pj) 

1023  kg/m3 

1000  kg/m3 

968  kg/m3 

Viscosity  of  pore  fluid  (i/) 

0.00105  kg/m-s 

0.001  kg/ms 

0.096  kg/m  s 

Porosity  (/?) 

0.385 

0.384 

0.36 

Permeability  (k) 

2.5  xlO-11  m2 

1.1  x  10"10  m2 

8.43  xlO"11  m2 

Tortuosity  (a) 

1.35 

1.45 

1.47 

Material  exponent  (n  =  m) 

0.0901 

0.0467 

0.3902 

Rigidity  coefficient  (77*  =  7 p  +  (4/3)75) 

2.6823  xlO8  Pa 

3.871  xlO8 

1.472  xlO6  Pa 
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TABLE  II:  Material  exponents  and  compressional  rigidity  coefficients  used  in  the  Biot- 
Buckingham  model  for  each  sediment.  The  parameters  were  determined  from  a  fit  to  the  data 
allowing  both  parameters  to  be  free  in  the  first  case  and  fixing  n  =  1  in  the  second  case. 


Sediment 

Fit  with  two  free  parameters 

Fit  with  one 

free  parameter 

n 

7 v 

n 

7 v 

SAX99 

1.088 

41.96  Pa 

1.00 

146.82  Pa 

Water/Glass  Beads 

0.915 

182.94  Pa 

1.00 

54.79  Pa 

SilOO/Glass  Beads 

1.059 

57.89  Pa 

1.00 

115.55  Pa 
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FIG.  1:  (a)  Sound  speed  and  (b)  attenuation  measured  during  SAX99.  Details  of  these  measure¬ 
ments  can  be  found  in  Ref.  12,  Also  shown  are  the  sound  speed  and  attenuation  predictions  for 
Biot  Theory  (solid  line)  and  the  Buckingham  model  (dashed  line)  using  the  parameters  given  in 
Table  I. 

FIG.  2:  Sound  speed  measurements  made  in  a  glass  bead  sediment  saturated  with  (a)  water 
and  (b)  silicone  oil.  Attenuation  measurements  made  in  a  glass  bead  sediment  saturated  with 
(c)  water  and  (d)  silicone  oil.  Also  shown  are  the  predictions  of  the  EDFM  (solid  line)  and  the 
Buckingham  model  (dashed  line)  using  the  parameters  given  in  Table  I.  The  error  bars  on  the 
data  points  represent  uncertainties  in  the  position  of  the  hydrophones  as  well  as  the  distribution 
of  multiple  measurements  in  each  sediment. 

FIG.  3:  (a)  Shear  wave  speed  and  (b)  attenuation.  The  data  point  at  1  kHz  is  a  shear  wave  mea¬ 
surement  made  during  SAX99  by  Richardson11.  The  shear  wave  speed  and  attenuation  predicted 
by  Biot  Theory  (dashed  line),  the  Buckingham  model  (dotted  line),  and  the  Biot-Buckingham 
model  (solid  line).  Each  of  these  models  was  fit  to  the  data  point. 


FIG.  4:  (a-c)  Sound  speed  and  (d-f)  attenuation  predictions  using  the  Biot-Buckingham  model 
for  the  sediment  at  SAX99,  the  water  and  glass  bead  sediment,  and  the  silicone  oil  and  glass  bead 
sediment.  The  grain-to-grain  shearing  parameters  in  the  Biot-Buckingham  model  were  determined 
in  two  ways.  For  the  solid  line,  the  values  of  both  the  compressional  rigidity  coefficient,  7P,  and 
the  material  exponent,  n,  which  gave  the  best  fit  to  both  the  sound  speed  and  attenuation  data 
were  found.  For  the  dashed  line,  the  material  exponent  was  set  to  n  =  1  while  the  value  of  was 
determined  from  the  best  fit. 
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Hefner  and  Williams,  Fig.  3 
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Hefner  and  Williams,  Fig.  4 


